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It has been shown recently that the Jarzynski equality is generalized under nonequilibrium feed- 
back control [T. Sagawa and M. Ueda, Phys. Rev. Lett. 104, 090602 (2010)]. The presence of 
feedback control in physical systems should modify both Jarzynski equality and detailed fluctuation 
theorem [K. H. Kim and H. Qian, Phys. Rev. E 75, 022102 (2007)]. However, the generalized 
Jarzynski equality under forward feedback control has been proved by consider that the physical 
systems under feedback control should locally satisfies the detailed fluctuation theorem. We use the 
same formalism and derive the generalized detailed fluctuation theorem under nonequilibrium feed- 
back control. It is well known that the exponential average in one direction limits the calculation of 
precise free energy differences. The knowledge of measurements from both directions usually gives 
improved results. In this aspect, the generalized detailed fluctuation theorem can be very useful in 
free energy calculations for system driven under nonequilibrium feedback control. 
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Recent development in statistical physics relationships 
describe nonequilibrium dynamics in terms of equali- 
ties [l|-|3|. In particular, the Jarzynski's equality [1| 
and the Crooks detailed fluctuation theorem [2| can be 
used to calculate equilibrium free energy differences from 
nonequilibrium work distribution. A system initially at 
equilibrium with temperature (inverse) /3 — 1/fceT (ks 
is the Boltzmann constant) is externally driven from its 
initial state to final state by nonequilibrum process sat- 
isfies the detailed fluctuation theorem P{W)/P(—W) = 
exp(/3(W — AF)) and its integrated version, the Jarzyn- 
ski equality (exp(— j3W)} = exp(— /3AF). Where W de- 
notes the work performed on the system, AF is the free 
energy difference of the system between its final and ini- 
tial equilibrium states and P(±W) is the work probabil- 
ity distribution in forward (+) and reverse (-) direction. 
The (exponential) average (.) is taken over the ensemble 
of nonequilibrium trajectories. These relationships has 
been verified in experiments [J, [5j as well as simulations 
[a, H and widely used in many branches of Science (see, 
eg.jg). 

Various experiments and simulations has been per- 
formed by adopting a suitable time-dependent driving 
scheme described by an external control switching proto- 
col. Even though the Jarzynski equality and the detailed 
fluctuation theorem are valid for any time-dependent 
driving scheme, the efficiency of a nonequilibrium switch- 
ing simulations which use these relationships to estimate 
precise free energy difference is depends on the switch- 
ing function [9( . The practical difficulties faced in precise 
calculation of free energy differences in simulations [101 ] 
initiated further developments in these theories 11] . In 



our recent work, we have discussed the analogy between 
the optimized switching free energy simulation and the 
nonequilibrium systems which are subjected to a feed- 
back control 12 1 . 

The evolution of the physical systems can be modified 
or controlled by repeated operation of an external agent 
called controller [LJ, [lj] . In contrast to open loop con- 
troller which operates on the system blindly, the feedback 
or closed loop controllers use information about the state 
of the system. The feedback is the process performed by 
the controller of measuring the system, deciding on the 
action given the measurement output, and acting on the 
system [15|. For example, in a single molecule Atomic 
Force Microscopy experiment, the external agent is an 
electric feedback circuit which detects the motion of the 
cantilever and manipulate the control force proportional 
to its velocity [16j j . The proper utilization of the informa- 
tion about the state of the system in feedback control ef- 
fectively improves the system performance 13l4l7lj . How- 
ever, the presence of feedback control in physical system 
modifies both the Jarzynski equality and the fluctuation 
theorem 



m- 



Recently, the Jarzynski equality is generalized to an 
experimental condition in which the system is driven be- 
tween two equilibrium state via nonequilibrium process 
under forward feedback control 18[ . Since the work that 
performed on a thermodynamic system can be lowered 
by feedback control [18[ , this feedback mechanism can be 
helpful in simulation for sampling rare trajectories and 
calculate precise free energy differences [L2| . The equilib- 
rium free energy difference for the driven system (which 
locally satisfies the detailed fluctuation theorem) under 
nonequilibrium feedback control in forward direction can 
be calculated from the generalized Jarzynski equality [181 ] 
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where a = (3(W — AF) and / is the mutual information 
measure obtained by the feedback controller [18]. The 
average is taken from the work distribution in forward 
direction with feedback control. 

The second law of thermodynamics can be quantita- 
tively described by the fluctuation theorem which are 
closely related to the Jarzynski equality. If the experi- 
ments or simulations has been performed under feedback 
control, both the Jarzynski e qual ity and the fluctuation 
theorem should be extended 16]. However, the gener- 
alized Jarzynski equality (EqfT]) under forward feedback 
control has been proved by consider that the physical sys- 
tems under feedback control should locally satisfies the 
detailed fluctuation theorem |18j. In this letter, we use 
the same formalism and derive the generalized detailed 
fluctuation theorem under nonequilibrium feedback con- 
trol. In order to calculate the free energy differences pre- 
cisely in simulations one require information of work dis- 
tribution in both forward and reverse direction 191-421 1 . 
In this aspect, one needs the generalized detailed fluctu- 
taion thoerem under nonequilibrium feedback control. 

The feedback control enhances our controllability of 
small thermodynamics systems 13H17]. At a given time, 
the controller measure the partial state of the system. 
The result of the measurement determines the action the 
control will take. The additional information on the sys- 
tem provided by the measurement further determines 
the system states. Suppose, the controller perform a 
measurements on a stochastic thermodynamics system at 
time t m . Let T m be the phase-space point of the system 
at that time, P[F m ] its probability and y the measure- 
ment. Depends on the controller measurements the mea- 
surement outcome y occurs with probability P[y] [18j. 
The information obtained by the controller measurement 
can be characterized by the mutual (feedback) informa- 
tion measure 



U 



I[y, r m ] = h 



PWm\ 

P[y] 



(2) 



where P[y|r m ] is the conditional probability of obtaining 
outcome y on condition that the state of the system is 
r m . The above equation is rewritten as, 



JlvX m ] _ 



P[yK 
P[y] 



(3) 



In experiments and simulations, the free energy differ- 
ence between the two equilibrium states can be calcu- 
lated in general by pulling the sytem from one equilib- 
rium state to another state along a switching path. The 
path connecting the two states in the time period r will 
be parameterized using the variable A, with < A < 1. 
The switching rate describes the nature of the switching 
process to be an equilibrium (infinitely slow) or nonequi- 
librium (fast). If the experiments performed under feed- 
back control, the switching control parameter A depends 



on the outcome y after t m 18]. That is, whenever the 
controller made measurements, there is a corresponding 
changes in the switching parameter for next time step, 
which is denoted as A( t;y ). Between the every stages of 
the controller measurements, the value of outcome y is 
fixed and the corresponding switching paremeter \(t;y) 
doesnot change. In this time interval for each stages of 
X(t-,y), we consider the system should locally satisfies the 
detailed fluctuation theorem 11811 . 
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where P\, t , y) [r(i)] is the probability of obtaining the 
outcome y in forward direction with switching protocol 
A( t .j,\ and P,t [T^t)] is the probability of obtaining the 

same outcome y [18J in reversed direction of phase point, 
r^(t), with corresponding switching protocol AL. ■>. Here, 
cr[r(i)] is the work value obtained in the forward direc- 
tion and its time reversal work value 



a[T\t)} = -cr[T(t)}. 



(5) 



Let P F [X] = P F [a[T],I[y, T]} be the the joint probabil- 
ity of obtaining the work value cr[T] for a given feedback 
information measure I[y, F] of measurement outcome y 
in forward direction. Due to the repeated measurements 
of the controller, the particular measurement outcome 
y may occurs in various sta ges of the experiment. In 
analogous to previous studies [23|,|2J], this (joint) proba- 
bility can be obtained from the nonequilibrium ensemble 
of variable (forward) switching trajectories as, 



p F [x] = J p[ y '\T' m ]p Xlt . yl) \r(t)W\y',r' m ]-mT]) 

5(a[T(t)}-a[f])dy'D[T(t)}, (6) 

where S(x) is the Dirac delta function which has a prop- 
erty 5(— x) = 8{x). It should be noted that y' and T' m in 
the conditional probabilities are dummy variables and it 
has appropriate values for each outcome of the controller 
measurements, see Eq.([3]). 

Combining Eq.([3]) and Eq.(j4]), then Eq.© becomes, 



P F [X] 
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5{I[y',T' m ]-I[y,T])8{a[T{t)\-a[T]) 
dy' D[T(t)), 



P F [X] = e°^ +I ^ / P[j/']P,t [Ft(i)] 



8(I[y',T' m ]-I[y,T})5(a[T(t)}-a[T}) 
dy' D[T(t)]. (7) 



Let Pr[X^] = P R [-a[T},I[tf,f]] be the the joint prob- 
ability of obtaining the work value — <r[r] for a given feed- 
back information measure 7[yT,r] of measurement out- 
come \p in reverse direction. This (joint) probability can 
be obtained from the nonequilibrium ensemble of variable 
(reverse) switching trajectories as, 



Pr[^] 



P[y'\r'jp x i . [rt(t)] 
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S(I[y'X m ]-I[^,f])6(a[rHt)}+a[T]) 

dy' Dp(t)}. (8) 

To this end we will use above equations and derive the 
generalized detailed fluctuation theorem under feedback 
control either in forward direction [18] or both directions 
as follows. 

If the system has same feedback control in both direc- 
tions then under the controller measurements condition 
Iffi , T] = I[y, T] and using Eq.©, we can rewrite Eq.© 
as, 

P R [X*] = e 7 ^ [p[y']P x i [rt(t)] 



6(I[y',T' m ]-I[y,T])6(a[T\t)] + a[T}) 

dy' D[T\t)}. (9) 

Since S(-x) = S(x) and D[T^{t)] = D[T(t)} 18], com- 
bining Eq.© and Eq.© in Eq.© we can obtain the 
generalized detailed fluctuation theorem under feedback 
control in both direction as 

P F [X] = e°ft +I VfiPR[Xi]e-Wfi, 



Pr[X1] 
The above equation can be written simply as 
Pf[ct,I] 



Pr[-v,I] 



(10) 



(11) 



Our result shows that for a given feedback information 
measure / in both directions, the physical system under 
feedback control satisfies the detailed fluctuation theo- 
— £1ZL — e a . Since the feedback control enhances 
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our controllability of small thermodynamics systems, the 
proper choice of feedback mechanism in free energy sim- 
ulations can be useful for precise free energy estimates 
instead looking for optimized switching protocols [Ijj . 

If the system has feedback control only in forward di- 
rection, the generalized detailed fluctuation theorem un- 
der forward feedback control can be obtained from the 
following reverse experimental condition. Based on the 
informations about switching protocols for each outcome 
in forward feedback control experiment, without feedback 



control, we can perform same variable switching proto- 
cols experiment in reverse direction 18]. This is equiva- 
lent as an open loop controller which operates on the sys- 
tem in reverse direction. In such a case, P[y'|r^] = P[y'] 
in reverse direction and the open loop controller implic- 
itly has information about the forward feedback infor- 
mation measure I[y, T] for corresponding switching pa- 
rameter, AL ,n, in reverse direction, see Eqs.© and ©. 
Then, we can write Eq. © as 
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6(I[y',T' m }-I[y,T})5(a[T\t)}+a[T}) 
dy' D[T\t)}. (12) 

As similar to earlier derivation, combining Eq.© and 
Ea. dl^]) in Eq.© we can obtain the generalized detailed 
fluctuation theorem under forward feedback control as 



Pf[X] 



Q a[T]+I[y,r] 



Pr\x^] 
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a a[r]+i[y,r] 



(13) 



In order to prove the generalized Jarzynski equality 
under forward feedback control, we measure the quantity, 

(e- ff - J ) = f P F [X] e - a F]-ily,r] d ~ ^ 



From Eq. ([13 
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f P R [X^]dydT 



(14) 



we obtained the generalized Jarzynski equality under for- 
ward feedback control [18j. In order to get more insight 
of forward mutual information measure, we calculate the 
quantity [18| 



(e- a ) 



P F [X] e- a[ii dy dT 



Using Eq. ()13[) , the average of above equation can be writ- 
ten as 

(e- a ) = f P R [X^] e^Uy df 



(15) 



= / P R [X* } dy dT, 



7, 



where 7 = J P R [X* ] dy dT is the feedback control char- 
acteristics which is a measure of the correlation between 
the dissipation and the information [18j and P R [X* ] — 



Pji[X^]e I ^' r ^ is the special case 18 1 of the joint proba- 
bility distribution in reverse direction. 

Even though there is no feedback control in reverse di- 
rection, due to the implementation of variable (reverse) 
switching protocols in accordance with forward feedback 
control experiment, the forward mutual information mea- 
sure can also obtained from the reverse direction as 



mr\ 



In 



p R [xn 



(16) 



In conclusion, we have derived the generalized detailed 
fluctuation theorem under nonequilibrium feedback con- 
trol. It is well known that the exponential average in 
one direction limits the accurate calculation of free en- 
ergy differences in simulation. The knowledge of mea- 
surements from both directions usually gives improved 
results. Thus, the generalized detailed fluctuation the- 
orem can be very useful in free energy simulation for 
system driven under nonequilibrium feedback control. 
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